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This  report  examines,  from  the  viewpoint  of  singular  perturbation  theory,  a  fourth 
order  partial  differential  equation  which  was  derived  by  Paidoussis  [  1  ]  as  a  model  for  the 
behavior  of  a  thin  cylinder  in  an  axial  flow.  It  is  found  that  for  sufficiently  large  form  drag 
(C{  <  Vx)  and  small  flexural  rigidity  the  influence  of  the  higher  order  boundary  conditions  is 
restricted  to  the  boundary;  i.e.,  the  reduced  equation  is  a  good  approximation.  Further¬ 
more  for  small  frequencies  the  downstream  boundary  conditions  can  be  ignored  in  the  sense 
that  outside  of  the  very  end  of  the  cylinder,  their  effect  on  the  solution  is  negligible. 

Finally,  an  examination  of  the  characteristics  of  the  reduced  PDE  leads  one  to  conjecture 
that  this  remains  true  in  the  case  cj  <  Vx. 
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I.  INTRODUCTION 

This  report  examines,  from  the  viewpoint  of  singular  perturbation  theory,  a  fourth 
order  partial  differential  equation  which  was  derived  by  Paidoussis  f  1  ]  as  a  model  for  the 
behavior  of  a  long  thin  flexible  cylinder  in  an  axial  flow.  The  motivation  for  our  study  of 
this  phenomenon  as  well  as  a  summary  of  previous  work  in  the  area  may  be  found  in  [2] . 

In  particular,  Ortloff  and  Ives  [3]  and  Kennedy  [2,4]  succeeded  in  finding  closed  form 
solutions  by  passing  to  the  limit  of  zero  flexural  rigidity,  thereby  reducing  the  governing 
equation  to  second  order.  We  shall  primarily  be  concerned  with  the  validity  of  the  reduced 
equation,  determining  suitable  boundary  conditions,  and  assessing  the  degree  of  their 
influence  on  the  solution. 

The  complete  partial  differential  equation  (PDE)  is  given  by  (the  notation  is  that  of 
of  [2]): 

EIyXXXx  +  (M  +  m)ytt  +  2MUytx 

+  (MU1 2 3 4)(1  -  2ct  (L  -  x)/dc  -  2c^)yxx 

+  (MU2)(2cn/7rdc  +  2ct/dc)yx 

+  (MU)(2cn/7rdc)yt  =  0  (1.1) 

where 

y(x,t) 
x 

t 

M 
m 
U 
L 

dc 
ct 


E 
I 

1.  M.  P.  Paidoussis,  “Dynamics  of  Flexible  Slender  Cylinders  in  Axial  Flow.”  Journal  of  Fluid  Mechanics. 
Vol26, 1966,  pp  717-751. 

2.  R.  Kennedy  and  E.  Strahan,  “A  Linear  Theory  of  Transverse  Cable  Dynamics  at  Low  Frequencies  ” 
Naval  Underwater  Systems  Center  TR  6463,  June  1981 . 

3.  C.  R.  Ortloff  and  J.  Ives,  “On  the  Dynamic  Motion  of  a  Thin  Flexible  Cylinder  in  a  Viscous  Stream,” 
Journal  of  Fluid  Mechanics.  Vol  38,  part  4, 1969,  pp  713-720. 

4.  R.  M.  Kennedy,  “Low  Frequency  Transverse  Motion  Response  of  Cable  Towed  Array  Systems,”  NUSC 
Technical  Memorandum  No.  93231,  Naval  Underwater  Systems  Center,  New  London,  CT,  July  1979. 


=  crosstrack  position  of  the  cylinder  measured  at  x 
=  along-track  space  variable  increasing  downstream 
=  time 

=  virtual  mass  of  fluid  per  unit,  7rp(dc/2)2,  where  p  is  the  fluid  density  [  1  ] 
=  mass  of  the  cylinder  per  unit  length 
=  free  stream  velocity 
=  cylinder  length 
=  cylinder  diameter 

=  tangential  (longitudinal)  drag  coefficient 

=  drag  coefficient  added  to  account  for  a  drogue  and/or  form  drag  at  the 
end  of  the  cable 

=  coefficient  of  the  linearized  normal  drag  expression 

=  Young’s  modulus 
=  moment  of  inertia  =  7r/4(dc/2)^. 
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TOW  POINT 


Figure  1 .  Physical  model  approximated  by  equation  (1.1). 

Several  approximations  were  involved  in  the  derivation  [  1  ]  of  equation  (1.1).  These 
include  a  hydrodynamic  approximation  which  uses  the  virtual  mass  M  =  irp( dc/2)2  concept 
(valid  for  a  straight  cylinder)  to  calculate  the  normal  component  of  the  fluid  momentum 
of  a  curved  cylinder,  and  a  mathematical  approximation  which  assumes  yx,  the  tangent  of 
the  angle  of  the  cylinder  with  the  x-axis,  to  be  sufficiently  small.  The  latter  is  reflected  in 
several  places:  the  neglect  of  nonlinear  drag  terms  in  the  first  order  derivatives,  the  term 
“1”  in  the  coefficient  of  yxx  which  approximates  the  cosine  of  the  angle  of  cylinder’s  axis, 
and  the  use  of  x  =  L  for  the  end  point  of  the  cylinder. 

A  dimensionless  version  of  ( 1 . 1 )  is  given  by 


where 


e3hl?&££  +  67Jtt  +  +  b(€  -  1  "  eq )Vtf  +  (b  +  d)rj£  +  drjT  =  0 

{  -  f, r=  ft, ,({,»=  • 

L  L  A 

A  =  forcing  function  amplitude 
e  =  dc/L 

0  =  M/(m  +  M) 

h  =  ±JL 

16  p  u2 

b  =  20ct 
q  =  (2c't-l)/2ct 

d  =  20c  n/tr 


d.2a] 


(1.2b; 

d.2c] 
( 1 -2d] 


(1.2h 


In  the  study  that  follows,  we  shall  be  interested  in  the  behavior  of  the  solution  as  e^h  -*■  0. 
One  way  to  achieve  this  is  to  let  the  relative  cylinder  diameter  approach  zero;  however,  it 
should  be  noted  that  0  depends  on  dc  through  M  and,  possibly,  through  m.  We  shall  assume 


that  the  mass  of  the  cylinder  per  unit  length  is  proportional  to  so  that  0  =  constant.  This 
includes  the  primary  case  of  interest,  that  of  neutral  buoyancy  for  which  m  =  pir(dc/2)^,  and 
thus  0  =  Vi. 

The  boundary  conditions  of  (1.2)  are  those  of  a  cylinder  pinned  at  one  end  (£  =  0) 
and  free  at  the  other  end  (£  =  1).  The  boundary  conditions  at  the  pinned  end  are  provided 


by  a  forcing  function  <f>  j 

t}(0,t)  =  0j(r)  (1.3) 

and  a  requirement  that  the  bending  moment  be  zero  [  1  ] 

T?££(0,T)  =  0  .  (1.4) 

Similarly,  there  must  be  no  bending  moment  at  the  free  end 

t?^(1,t)  =  0  .  (1.5) 

The  second  boundary  condition  represents  a  balancing  of  the  lateral  shear  force  and  is  given 
by  the  somewhat  complex  expression  [  1  ] 

El  y xxx  +  fMU(y t  +  Uyx)  -  (m  +  fM)xeytt  =  0  (1.6) 

where  f  is  a  constant  less  than  1  and 
L 

Xe=A  /  A(x)dx  0-7) 


with  A  the  cross-sectioned  area  for  an  end  tapered  from  L  -  £  to  L. 

To  simplify  our  computations  we  shall  take  f  =  1  and  assume  that  £«  L  so  that 
xe  =  0.  However,  note  that  in  the  case  of  a  drogue,  which  in  [2]  is  modeled  as  an  extra 
piece  of  cylinder,  xe  is  not  zero.  Alternatively,  a  drogue  may  be  considered  a  point  mass 
mg  so  that  the  last  term  is  replaced  by  the  acceleration  mpytt.  In  dimensionless  units 
this  gives 

he2r?£££  +  0(j?t  +  r?£)  +  o0i7TT  =  0  (1.8) 

where  a  =  mass  of  drogue/mass  of  cylinder  at  neutral  buoyancy.  (More  precisely, 
a  =  mD/ML.) 

As  the  flexural  rigidity  goes  to  zero,  (1.8)  becomes  an  expression  of  the  first  order 
in  f .  As  will  be  detailed  in  section  II,  this  implies  that  the  boundary  condition 

0?T  +  »?$)  +  ocnTT  =  0  ,  (1.9) 

derived  from  (1.8),  is  appropriate  for  the  reduced  equation.  Thus,  one  must  resist  the  tendency 
to  dismiss  or  remain  indifferent  to  (1.6)  simply  because  it  is  the  higher  order  boundary 
condition. 

Since  the  problem  as  stated  here  is  linear,  we  shall  find  it  expedient  in  a  great  portion 
of  the  development  to  use  separation  of  variables.  For  this  purpose  tj  is  allowed  to  assume 
complex  values,  and  it  is  understood  that  the  physical  solution  is  given  by  its  real  part.  The 
driving  force  is  taken  as  0  j  (r)  =  e1TCJ  so  that  real  (0j )  =  cos  (wr),  and  more  general  solutions 
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may  be  obtained  by  super-position.  We  use  this  approach  in  section  II  and  solve  the 
resulting  ordinary  differential  equation  (ODE)  by  singular  perturbation  analysis.  Many  of 
the  specific  computations  rest  on  the  assumption  oo  «  1 ;  i.e.,  are  for  large  wavelengths. 
However,  this  does  not  preclude  (and  in  fact  would  even  indicate)  the  qualitative  validity 
of  our  results  in  a  more  general  context.  The  mathematically  crucial  assumption  that 
c't  >  Vi  is  abandoned  in  section  III,  but  we  are  still  able  to  derive  several  interesting  although 
incomplete  results.  In  particular,  the  reduced  PDE  (as  opposed  to  ODE  of  section  II)  is 
examined  by  studying  its  characteristics.  Only  the  “steady  state’’  solution  is  considered; 
aspects  relating  to  initial  conditions  and/or  stability  are  conveniently  ignored.  Our  method¬ 
ology  shall  not  be  uniformly  rigorous  but  is  intended  to  give  a  reasonably  logical  analysis 
while  pointing  out  the  various  gaps  that  remain  to  be  filled. 


II.  A  SINGULAR  PERTURBATION  ANALYSIS  FOR  THE  CASE  c[  >  K 

This  section  is  devoted  to  a  singular  perturbation  analysis  in  which  the  effects  of  the 
structural  rigidity,  as  reflected  in  the  coefficient  e^h  of  eq.  (1.2a),  are  small.  Following 
O’Malley  [5] ,  we  define  a  perturbation  parameter 

H=y/e*h  (2.1) 


and  try  to  find  an  asymptotic  expansion  in  that  parameter.  Note  that  equation  (1.2a)  really 
involves  two  small  parameters,  p  and  e,  corresponding  to  the  structural  rigidity  and  relative 
cylinder  diameter,  respectively.  This  aspect  shall  be  addressed  in  more  detail  later. 

We  separate  variables  making  the  substitution 

t ?(€,T)  =  eiwTv«)  (2.2) 

and  obtain  the  fourth  order  ordinary  differential  equation  (cf.,  (A.l),  (A.4),  and  (A.5)) 

p2v( 4)  +  a-j(£  _  l  _  eq)  v^)  +  -hc^v  =  0  (2.3a) 


where  0  <  £  <  1  and 
o']  =  b 

<*2  =  b  +  d  +  20eico  -  b(  1  -  v) 

4oj  [ id  eco\  A  2 

“3=TVi  ■  T/-bo  • 

The  notation  v^1^  is  used  to  represent  the  j1*1  derivative  of  v  with  respect  to  £. 
Following  [5] ,  we  seek  a  composite  asymptotic  expansion  of  the  form 


(2.3b) 


5.  R.  O’Malley,  Jr..  Introduction  to  Singular  Perturbations,  Academic  Press.  NY.  1974,  Chapters  3 
and  8. 
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v($,p)  =  Gl($)  +  G2(S)  +  M[Gu(S)  +  G21(*)]  +  ... 

+  p2[L0(*)  +  pL1($)  +  ...]  e'gl(5)/M 

+  p2[R0«)  +  pR1a)  +  ...]  e"g2(*)/M 

(2.4) 

with  boundary  conditions 

v(0,p)  =  1 

(2.5a) 

O 

II 

/— s 

=t 

O" 

Ml* 

>* 

(2.5b) 

p2v(3)(  j  ^  +  0e(V£(  1  ,p)  +  icjv(  1  ,p))  =  0 

(2.5c) 

v^(l,p)  =  0 

(2.5d) 

The  first  line  of  equation  (2.4)  contains  the  so-called  reduced  solution  G|(£)  +  G2(£)  of 
(2.2),  obtained  by  setting  p  =  0.  Since  the  order  of  the  resulting  equation  is  2,  there  are  two 
linearly  independent  solutions,  Gj  and  G2,  determined  by  the  boundary  conditions  (2.5a) 
and  (2.5c).  The  functions  Gy  are  the  higher  order  terms  in  the  expansion  modulo  the 

“boundary  layer”  terms  found  in  the  next  two  lines.  For  small  p  the  exponential  terms 
-g  |  fax  -g^  fax 

e  1  and  e  z  are  significant  only  in  a  small  area,  so-called  inner  regions,  about  the  left 
hand  and  right  hand  boundaries,  respectively.  The  factors  of  p2  preceding  these  terms  are 
dictated  by  the  boundary  conditions  (2.5b)  and  (2.5d)  in  which  the  second  derivatives  of 
v  produce  a  factor  (1/p)2  which  must  be  canceled  by  the  coefficient  p2. 

The  boundary  conditions  (2.5a)- (2. 5d)  were  obtained  from  (1.3),  (1.4),  (1.5),  and 
(1.6)  -  (1.8)  through  the  separation  of  variables  (2.2).  Note  that  the  form  of  (2.2),  which 
in  conjunction  with  (2.5a)  also  implies  the  boundary  condition  77(0,7)  =  eitJT,  represents  no 
real  restriction  since  equation  (1.2a)  is  linear  in  rj  and  we  may  apply  superposition  to  con¬ 
struct  a  more  general  solution.  For  computational  simplicity,  we  have  set  a  =  0  in  equation 
(1.8).  Actually,  this  should  be  a  valid  approximation  (cf .  section  I)  in  the  absence  of  a 
drogue,  with  a  drogue  of  sufficiently  small  mass,  or  for  very  small  to. 

The  unknown  functions  of  (2.4)  are  determined  by  substituting  (2.4)  into  (2.3) 

and  equating  the  coefficients  of  powers  of  p  (or  of  pJe’^  since  approaches  zero 
faster  than  any  power  of  p)  to  zero.  This  leads  to  a  sequence  of  differential  equations  whose 
constants  are  determined  by  the  boundary  conditions  (2.5).  At  stage  j  a  set  of  simultaneous 
equations  (2.5a)  and  (2.5c)  specify  G jj  and  G2j.  Lj  and  Rj  are  then  determined  from  (2.5b) 
and  (2.5d)  respectively.  It  should  be  noted  that  the  coefficient  of  pi  in  (2.5c)  only  involves 
R^  for  k  <  j  since  the  coefficient  of  Rj  involves  powers  of  p  of  the  order  2  +  (2  +  j  -  3) 

=  j  +  1  and  higher.  We  begin  our  computations  with  the  determination  of  G  j  and  G2,  the 
reduced  problem. 

REDUCED  PROBLEM 

The  reduced  problem,  obtained  by  setting  p  =  0  in  (2.3)  is  given  by 

aj($-  1  -eq)v^  +  a2^^  +  a3^  =  0  (2.6a) 


with 

v(0)  =  1 


(2.6b) 
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v€(  1 )  +  iu/v(  1 )  =  0. 


(2.6c) 


This  statement  is  somewhat  unusual  in  that  the  boundary  condition  (2.6c)  is  derived  from 
the  highest  order  condition  of  the  full  equation.  (It  occurs  because  (2.5c)  involves  p,  and 
its  reduced  order  is  less  than  that  of  (2.5d).)  Nevertheless,  such  behavior  does  not  appear 
to  create  any  problems,  and  suitable  generalizations  along  the  lines  of  [6] -[7]  should  lead 
to  the  standard  theorem  [  5  ] 

lim  v(jLt,f)  =  v(£)  in  (0,1)  (2.7) 

p  0 

provided  the  coefficient  of  v^)  in  (2.6a)  is  uniformly  bounded  away  from  zero;  i.e..  orovided 

eq  >  6  >  0.  (2.8) 

It  might  seem  that  q  >  0  would  constitute  a  sufficient  condition  for  (2.8)  to  hold,  but  some 
extra  care  is  warranted.  Recall  that  p  =  (he2)1/2.  Strictly  speaking,  the  limit  in  (2.7),  as 
described  in  [5] ,  assumes  that  e  is  a  constant  since  it  also  appears  in  several  coefficients  of 
(2.3):  we  have  included  e  in  (2.8)  to  emphasize  this  point.  One  would  strongly  suspect, 
however,  that  (2.7)  remains  valid  as  long  as  (2.8)  is  true  even  if  e  ->  0  (and  thus  q-*») 
although  a  rigorous  demonstration  would  require  a  careful  examination  of  the  proofs 
appearing  in  [5]  -  [7] . 

A  closed  form  solution  to  (2.6a)  was  derived  in  [3)  in  terms  of  Bessel  functions  [8] , 
which  we  now  write  in  the  following  form  (cf.  Appendix,  eqs.  (A.9)-(A.  10)): 


Gi  +  G2  =  AF  j(b^z2)  +  BF2(b^z2) 


(2.9) 


where  A  and  B  are  constants, 

o°  (_i  )m 


3<b0z2 


'■2 


(¥)' 


m=0 


m!r(-p  +  m  +  1) 


(2.10) 


F2(b5z2)  = 


(-r 

mTfp  +  m  +  1) 

m=0 


(2.11) 


z“  =  b($  -  1  -  eq) 


(2.12) 


and  v  and  bQ  are  given  by  (A.4)  and  (A. 5).  The  reduced  boundary  conditions  (2.6b)  and 
(2.6c)  become 


AF|(bgz2)  +  BF2(b^)  =  1 


Ab^bF'|(bgz|)  +  BbgbF'^(bgzy)  +  icu[  AFj(bgzy)  +  BF2(bgzy))  =  0 


(2.13) 

(2.14) 


6.  G.  Birkhoff,  “On  the  Asymptotic  Character  of  the  Solutions  of  Certain  Linear  Differential  Equations,” 
Trans.  AMS.  Vol  9,  April  1908,  pp  219-230. 

7.  Turritin,  H.,  “Asymptotic  Solutions  of  Certain  Ordinary  Differential  Equations  Associated  with  Multi¬ 
ple  Roots  of  the  Characteristic  Equation”  Amer.  J.  Math..  Vol  58.  1936.  pp  364-378. 

8.  G.  N.  Watson.  A  Treatise  on  the  Theory  of  Bessel  Functions.  Second  Edition.  Cambridge  University 
Press,  1966,  p  97. 
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with 


:§  =  -b(l  +eq) 

(2.15a) 

=  -beq 

(2.15b) 

and  F'  indicates  the  derivative  of  the  function  F  with  respect  to  its  argument. 

Since  T(r  +  1 )  =  rr(r),  the  ratios  of  the  second  terms  to  the  first  terms  in  the  series 
expansions  of  Fj  and  Fj  (eqs.  (A.9)-(A.12»  are  bounded  by  (real  part  of  v  negative) 


b2z2 

b0z 


II 


(2.16) 


Thus,  for  sufficiently  small  Ibbgeql  (and  v  not  too  close  to  an  integer),  we  may  approximate 
(2. 14)  by  substituting  the  first  terms  of  the  power  series  of  Fj  and  Fj 


'b0b 


ice 


AV{-v  +  2) 

+  B 


r(-y  + 1) 


-bob" 


4r>  + 1 )  \-?bboeq 


\  v~x  ,  icc 

y  n>  + 


i) 


--bb^eq 


■) 


=  0 


(2.17) 


Some  typical  values  of  the  parameters  involved  are  [2]  ct  =  .01,  c_  =  .2,  e  =  10-\  and 

__  _  .  .  .  . .  ....  ...  _ .  . 


c[  <  200.  In  that  case  we  see  that  for  u)  «  10^  (from  eqs.  (A.4)  and  (A.5)) 


bb 


In - 4ioci>  . 


(2.18) 


Also,  expression  (2.16),  whosfe  value  is  approximately  |coeq|  for  z  =  Zj  and  |co|  for  z  =  zg,  is 
small  for 


oc  «  1 

Substituting  (2. 1 8)  in  (2. 1 7)  yields 
*'v  1 


+  B 


- - -  (iccyeq)1^  ( v~  +  icci'eq) 

ro>+n 


|T(-  is  +  2)  n-v+\) 

and  after  some  manipulation  the  approximation  +  icoveq  ~  gives 

B - -2r^-±J)  .  {iu)V€q)^  . 

A  iz2(-iz  +  i )  r(-*z  +  i  ) 


=  o 


(2.19) 


(2.20) 


(2.21) 


For  the  above  values,  as  long  as  v  is  not  close  to  a  negative  integer,  |B|  «  |A|.  How¬ 
ever,  the  physically  significant  quantity  is  not  the  relative  size  of  the  coefficients,  but  rather 
that  of  the  functions  Gj  =  AF  j  and  G2  =  BF2.  Taking  the  ratios  of  these  two  functions,  we 
find  from  (2. 10),  (2. 1 1),  and  (2.21)  that 
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G2(z) 

b  (V  )  r(-r>+  i) 

2 

t  w 

h2  2 
b0z 

G  j(z) 

A  \  4  /  r(v+  1) 

i>2(l  -v) 

uveq 

4 

2  <  I 

<  - - 1  \u>veq 

v2(l-v)\  I 

Since  real  ( v )  <  0  and  \z\2  is  a  minimum  for  z 2  =  z2.  Thus  we  find  that  for  oc  «  1 


(2.22) 


IG2($)I  «  iGjCf)!  Jin  [0,1]  .  (2.23) 

In  fact,  for  the  above  values,  v  ~  -7  and  eq  ~  .2  so  that 

IG2(J)I  <  .Olcc  IGj(J)I  .  (2  24) 


Intuitively,  what  has  occurred  is  that  F2  is  more  singular  than  F2  at  z  =  zj ;  however, 
Fj  and  Fj  are  about  the  same  magnitude.  The  boundary  condition  (2.14)  forces  BF2  to  be 
of  the  same  order  as  AFj ;  and  thus,  G2  =  BF2  is  much  smaller  than  G  j  =  AFj.  When  it  is 
appropriate  to  neglect  G->  (e.g.,  when  cj  «  1),  we  may  write  the  two  term  approximation 
of (2.13)  as 


1 _  (bQzo)  1 

T(-v  +  1 )  \  4  /  r(-t>  +  1 ) 

- 


(2.25) 


and  the  one  term  (zeroth  order)  as 


A  ~  f(-i>  +  1 ) 


(2.26) 


THE  INNER  REGION  EXPANSIONS 


In  this  section  we  seek  expressions  for  the  functions  g  j ,  g2,  Lj,  and  Rj  of  equation 
(2.4).  They  are  determined  by  substituting  (2.4)  into  equation  (2.3)  and  matching  powers 
of  p.  A  straightforward  computation  (Appendix  A)  reveals  that  the  most  singular  term, 
involving  vanishes  provided 

g  =  ±-j  -v/oj-  (-x)^/2  +  constant  (2.27) 

where  x  =  J  -1  -  eq.  In  order  that  our  expressions  remain  bounded  in  [0,1  ]  and  still  produce 
a  finite  contribution  at  the  boundaries  as  p  -*■  0  we  must  have  g j  (J  =  0)  =  g2(J  =  1 )  =  0,  and 
hence 

gl(l>  =  gl(*(D>«yvGT  [(1  +  eq)3/'2  -  (-x)3/2]  t2. 28a) 

g2(£)si2(x^  )=  1  [(_x)3/2  ■  (ec^)3/2]  •  (2-28b) 


Note  that  we  have  taken  the  liberty  of  using  the  same  symbol,  g.  to  represent  two  different 
functions,  g(J)  and  g(x).  The  dependency  g(x)  is  assumed  unless  otherwise  indicated  by  the 
context.  Also  x  is  not  the  same  as  the  variable  “x”  which  appears  in  equations  (1.1)  and 


Substituting  these  expressions  into  (2.3)  and  equating  the  coefficient  of  /z-iie_6/#i  to 
zero  yields  (eq.  (A-24)) 


Ro+s(l-^)Ro  =  0 


(2.29) 


which  has  the  solution 


I  (°H .  i\ 

t0  =  r0(-x)2  \al  2/ 


(2.30) 


with  exactly  the  same  equations  holding  for  Lg.  The  constant  rg  is  determined  by  satisfying 
the  boundary  condition  (2.5d)  up  to  order  0(1)  in  n  for  the  expansion  (2.4).  Making  use  of 

equation  (A.  1 6)  for  the  second  derivative  of  Rge  L  ,  we  have  up  to  0(  1 )  in  p 

G'j(l)  +  G5(l)  +  r0(g'2(-<q))2R0(-€q)  =  0  .  (2.31) 

(Recall  G  is  a  function  of  £,  whereas  g2  is  a  function  of  x.) 

2 

Once  again  making  the  assumption  that  Ibbgeql  is  sufficiently  small,  we  may  obtain 
an  approximate  solution  of  (2.3 1)  for  rg.  From  (2. 1 8),  (2.2 1 ),  (A.  1  lb),  and  (A.  1 2b)  we 
find  that 


G'jO) 

G  j  ( 1 ) 


--  -2I>+1)  (i-q)l- 

vz(-v+  i)r(-^+ 1)  hi 


2  r(-i>  +  3)  ,_i 

—  .  . :  (iw^eq) 

v  r(-^  +  i  ) 
2v(uJi'eqr^ 


and,  recalling  that  G'{  =  (bgb)2F'j\ 


(2.32) 


G^l)- 


r(-K  +  3)  eq 


(2.33) 


Thus,  we  may  ignore  G'j(l)  in  (2.31 )  and  approximate  G'2(  I )  by  (2.33).  Also,  since  =  b 
g2(-eq)  =  -y/b~ v/eq  n 


and 


l 

Ro(~e(i)  =  ro^eci)2 


Hence, 

2Av 2  ico  + 

T(-»  +  3)  eq 


r0(beq)(eq)  2 


(• 


0 


or 


2Av2 
r(-t>  +  3) 


(2.35) 

(2.36) 


(2.37) 


M2Roe 


As  in  the  case  of  the  reduced  equation,  we  are  interested  in  the  contribution  of 

-go  //LX 

relative  to  the  entire  solution.  Hence,  we  examine  the  ratio 


^Roe'82'" 

2  -  -  - 

2vL  r(-v+l)o)f  .2  4 

G1 

rx-,  +  3»  b  (eq) 

v_  3_ 

<-x>*  2  '  4  e-82'" 


v_  _  5 

-v  3 

—■  cj(eq)2  4 

D 

(-X)  2  "  4  e-82/« 

(2.38) 


Inasmuch  as  the  perturbation  parameter  p  is  assumed  small,  the  ratio  (2.38)  will  be  significant 
outside  the  right  hand  “boundary  region”  (i.e.,  a  neighborhood  of  x  =  -eq)  only  if  eq  is  very 
small.  Furthermore,  since  p  the  exponential  in  (2.38)  dominates  as  e  goes  to  zero, 

and 


lim  [ratio  (2.38)]  =  0  forx<-eq  .  (2.39) 

e  -*■  0 


Thus,  q  itself  will  have  to  be  very  small  if  (2.38)  is  to  be  0(1)  outside  the  boundary  region.  To 
be  precise,  let  us,  for  example,  set  the  edge  of  the  boundary  region  to  be  at  x  =  -eq  -  0. 1  and 
determine  the  magnitude  of  q  necessary  for  (2.38)  to  die  out  outside  the  region.  As  earlier, 
we  take  nominal  values  of  v  ~  -7,  e  =  10“^,  and  b  =  .01.  For  small  eq,  g2(-eq  -  .1) 
~\/b'(.l)^2  ~  2  •  10”^.  Then  (2.38)  will  be  0(1)  if  eq  «  1  and 

-19  J_l_ 

200w(eq)  4  (.1)  4  p2  e~002^  ~  1  .  (2.40) 

If  Young’s  Modulus  E  is  about  10^,  u  =  500,  and  p  =  1  in  cgs  units  [2] ,  then  h  =  0E/(16pU2) 
~  1/8  •  104  so  that  p  ~  3  •  10~^/2.  Substituting  these  values  into  (2.40),  we  obtain 
(since  u>  «  1 ) 


q< 


_J _ 

j0500  ’ 


(2.41) 
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i.e.,  q  has  to  be  essentially  zero  for  the  “inner”  solution  t  to  have  a  significant  effect  (except 
at  the  very  boundary). 

REMARKS  ON  CONVERGENCE 

For  eq  =  0,  the  coefficient  tq  of  (2.37)  becomes  infinite.  It  is  thus  natural  to  inquire 
as  to  conditions  (i.e.,  how  small  an  h  is  necessary  for  a  given  e  and  q)  such  that  the  first  term 
of  the  asymptotic  expansion  is  a  good  approximation  to  the  true  solution.  Such  a  study 
should  be  straightforward,  although  we  shall  not  undertake  it  in  the  present  paper. 

An  interesting,  but  more  difficult  question,  is  the  validity  of  (2.4)  as  an  asymptotic 
expansion  in  e.  Some  preliminary  computations  lead  us  to  conjecture  that  the  expansion  is 
not  uniform  in  e  as  e  -*■  0  (with  h  fixed)  on  the  full  interval  £e[0,l  ] ,  but  it  is  uniform  on  any 
closed  proper  subinterval  [0,1  -  5]  where  0  <  5.  (There  should  be  no  problem  at  the  left 
hand  boundary.)  One  might  also  consider  various  combinations  of  h,  e,  and  q  simultaneously 
approaching  zero  and  determine  conditions  on  their  relative  magnitudes  which  assure  the 
asymptotic  convergence  of  expansion  (2.4). 

CONCLUSIONS 

For  q  >  0  (i.e.,  c^  >  1  /2)  and  small  flexural  rigidity  (h  -►  0)  the  influence  of  the  higher 
order  boundary  conditions  is  restricted  to  the  boundary.  In  other  words,  the  reduced  equa¬ 
tion  is  a  good  approximation.  This  seems  to  imply  that  the  appropriate  downstream  boundary 
condition  in  the  absence  of  flexural  rigidity  is  ijT  +  r?f  lf=i  =  0  or,  more  generally,  (1.9).  This 
result  does  not  appear  obvious  from  direct  physical  considerations.  Furthermore,  for  small 
frequencies  (condition  2.19),  we  can  ignore  all  right  hand  boundary  conditions  in  the  sense 
that,  outside  of  the  very  end  of  the  cylinder,  the  effect  of  the  reduced  boundary  condition 
(1.9)  on  the  solution  is  negligible  (equation  (2.23)).  The  reader  is  warned,  however,  that 
these  comments  are  all  qualitative.  For  a  more  precise  statement,  he/she  should  examine  the 
mathematical  details  which  have  been  presented  throughout  this  section. 

III.  THE  CASE  c'  <1/2  (A  POTPOURRI  OF  CONJECTURES) 

The  analysis  of  the  previous  section  breaks  down  when  c't  <  1 12  (i.e.,  q  <  0).  In 
that  case  the  ODE  (2.3)  obtained  by  separation  of  variables  possesses  a  turning  point  at 
£  =  1  +  eq  so  that  the  composite  expansion  (2.4)  is  no  longer  valid,  and  the  possibility  of 
extremely  anomalous  behavior  cannot  be  ruled  out  ([5] ,  Chapter  8).  Furthermore,  an 
examination  on  the  reduced  PDE  (obtained  by  setting  h  =  0  in  (1.2a))  reveals  that  it  changes 
from  hyperbolic  type  to  elliptic  at  $  =  02e/b  +  1  +  eq;  hence,  an  elliptic  region  will  exist  if 
q  <  -02/b.  Of  major  concern,  then,  is  a  consistent  formulation  of  the  problem:  What  are 
the  appropriate  boundary  conditions  for  the  PDE  ( 1 .2a)?  Is  the  second  order  reduced 
equation  a  valid  approximation  to  (1.2a)?  and  What  are  suitable  boundary  conditions  for 
the  reduced  PDE?  A  detailed  investigation  of  these  issues  is  beyond  the  scope  of  the  present 
study.  Thus  we  shall  try  to  briefly  analyze  the  situation,  point  out  some  of  the  more  salient 
features,  and  close  with  a  list  of  conjectures  and  open  questions. 

Let  us  begin  with  an  examination  of  the  reduced  PDE 

^We  use  the  term  “inner  solution”  to  denote  a»2(Rq  +  + . .  >  although  strictly  speaking  we  are 

referring  to  the  part  of  the  composite  expansion  associated  with  the  inner  solution. 


ei?TT  +  20€»?j£  +  b(£  -  1  -  eq)i7^  +  (b  +  d)i?^  +  dr;T  =  0 


(3.1) 


which,  for  convenience,  we  rewrite  as 

t?TT  +  20t?t£  +  byq^  +  (b  +  d  ■  +  j-r?T  =  0  (3.2a) 

where  the  coefficient  y  is  a  function  of  £ 

7($)  =  -  q  •  (3.2b) 

This  equation  is  hyperbolic  for/J^-b7>0and  elliptic  for/3^-b7<0.  Since  £e[0,l  ] ,  a  neces¬ 
sary  condition  for  the  existence  of  an  elliptic  region  is  0  <  ]  (cf  equations  ( 1 ,2g)  and  ( 1 ,2h)) . 
The  width  of  such  an  elliptic  region  is  equal  to  e(  1  -  0  -  2c[)/2ct;  i.e.,  it  is  of  the  order  e  Thus, 
for  typical  values  of  e  ([21 ),  almost  the  entire  cylinder  will  lie  in  the  hyperbolic  region. 

CHARACTERISTICS  IN  THE  HYPERBOLIC  REGION 

The  dependency  of  the  solution  of  the  almost  linear  hyperbolic  equation  (3.2)  on 
the  initial  conditions  (r  =  0)  and  the  boundary  conditions  (where  £  =  0  or  1),  is  determined 
by  the  geometry  of  its  characteristics.  These  form  a  two  family  net  of  curves  in  the  £,  r  plane 
with  slopes  Xj(£)  and  X2<£)  satisfying  ([91,  [  10])  the  quadratic  equation 


b7X2  -  20X  +  1  =  0  . 

(3.3) 

That  is. 

*  _  0-v/02~b7 

X'  '  bX 

(3.4a) 

-X  _  0  +Vp2  -  by 

X2  ‘  b7 

(3.4b) 

Note  that  lim  Xj  =  1/20  and  lim  IX2I  =  00 . 
y  -*  0  7  -*•  0 


A  sketch  (not  exact,  but  qualitatively  correct)  of  these  characteristics  appears  in 
Figure  2.  Since  the  domain  of  dependence  of  a  given  point  is  determined  by  the  region 
bounded  by  the  characteristics  emanating  from  that  point  (see  Figure  3),  the  appropriate 
boundary  (initial)  conditions  for  equation  (3.2)  lie  on  the  axes  £  =  0,  t  =  0.^  More  precisely. 


we  must  specify  three  functions 

01  O')  =  t?(o,t) 

(3.5a) 

*|«)“»l(tO) 

(3.5b) 

02«>  =  t?£(°’$) 

(3.5c) 

t Recall  that  in  this  section,  we  make  the  assumption  c't  <  1/2  so  that  the  end  of  the  cable,  £  =  1 ,  lies  to  the 
right  of  7  =  0.  In  the  contrary  case,  we  have  the  situation  pictured  in  Figure  4,  and  the  domain  of  depen¬ 
dence  includes  the  boundary  £  =  1 ;  i.e.,  the  end  point.  If  this  were  not  true,  our  computations  would  be 
inconsistent  with  the  results  of  section  2. 
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Figure  4.  Example  showing  that  the  domain  of  dependence 
of  an  interior  point  P  includes  the  downstream  end  of  the 
cable  if  c't  >  'A  (i.e.,  7(1)  <  0). 

This  result  has  several  very  interesting  physical  implications.  If  the  solution  is  stable 
(i.e.,  if  the  transients  due  to  the  initial  conditions  4>  j  and  ^2  die  out),  then  the  steady  state  solu¬ 
tion  depends  only  on  the  tow  point  forcing  function  p  j .  This  appears  consistent  with  the 
results  obtained  in  section  II  for  c't  >  Vi; namely,  except  at  the  right  boundary,  the  effect  of  the 
downstream  boundary  condition  is  negligible.  Also,  the  vertical  characteristic  at  7  =  0 
represents  a  type  of  barrier.  Past  that  point,  disturbances  only  propagate  downstream.  The 
upstream  part  of  the  cylinder  (7  <  0)  does  not  “know”  what  the  downstream  end  (7  <  0)  is 
doing.  Finally,  at  7  =  0,  we  note  from  equations  (3.2b)  and  (1.1)  that 

MU2  =  2ct  iLu2(L-x)  ;  (3.6) 

3x2  dc  9x2 

i.e.,  the  transverse  force  of  the  fluid  momentum  exactly  balances  the  transverse  component 
of  the  tension  in  the  cylinder.  ^  At  such  a  point  the  cylinder  should  offer  vanishing  impedance 
to  an  external  vertical  force. 

On  the  other  hand,  these  observations,  in  turn,  generate  their  own  questions.  What  is 
so  special  about  the  drag  coefficient  value  cj  =  Vi,  and  does  not  the  zero  impedance  at  the 
barrier  imply  the  possible  occurrence  or  persistence  of  large  displacements,  i.e.,  instability? 

REMARKS  ON  THE  SMALL  ANGLE  APPROXIMATION 

We  propose  the  following  answer:  If  one  returns  to  the  original  derivation  of 
equation  ( 1 . 1 )  [  1  ] ,  it  is  found  that  q  should  actually  be  given  by 


1  /  _ l 

2Ct  \  y/\T 


+ Noted  by  R.  Kennedy,  personal  communication,  Naval  Underwater  Systems  Center. 
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but  that  a  small  angle  approximation,  1  A/1  ~  1 ,  was  made.  With  ( 1 ,2h),  y  vanishes 

in  [0,1  ]  if  and  only  if  c't  <  Vi.  On  the  other  hand  if  expression  (3.7)  is  used,  the  situation  is 
much  more  complex  since  the  condition  7(£)  =  0  may  depend  on  the  boundary  conditions 
through  the  term  containing  It  is  still  true  that  c[  <  Vi  represents  a  sufficient  condition 
for  7  to  vanish;  however,  the  “barrier”  position  now  becomes  a  function  of  time  (i.e.,  r).  It 
is  also  of  note  that  Figure  1  is  no  longer  a  valid  picture  of  the  characteristics.  We  cannot 
even  sketch  them  without  further  analysis  since  7  is  a  function  of  the  solution  through 
Nevertheless  it  seems  reasonable  to  assume  that  in  general,  there  will  be  no  vertical  character¬ 
istic,  and,  thus,  that  a  downstream  boundary  condition  must  be  supplied.  This  is  actually 
desirable  inasmuch  as  we  expect  the  full  fourth  order  equation  to  require  two  downstream 
boundary  conditions,  only  one  of  which  is  lost  in  going  to  the  reduced  PDE.  Furthermore, 
the  result  will  still  be  consistent  with  equation  (3.2)  and  its  lack  of  downstream  boundary 
conditions,  if  we  can  show  that  the  effects  of  the  nonlinearity  (3.7)  are  restricted  to  a  small 
region  near  the  boundary;  i.e.,  equation  (3.2)  and  its  solution  are  a  good  approximation. 

To  provide  a  convincing  argument  (and  quantitative  bounds)  one  would  have  to 
demonstrate  that,  given  any  boundary  and  initial  functions  (3.5a)  -  (3.5c)  such  that  1^1  is 
small,  then  the  slope  of  the  solution  of  (3.2),  ItjJ,  is  sufficiently  small  that  it  may  be  neglec¬ 
ted  in  (3.7)  except  possibly  near  the  boundary.  We  content  ourselves  with  a  heuristic 
argument  concerning  the  influence  of  the  downstream  boundary  condition  in  the  nonlinear 
(eq.  (3.7))  case. 

Suppose  that  for  all  r,  the  solution  r)  of  (3.1)  with  q  given  by  (3.7)  may  be  expanded 
in  a  Taylor  series  v^id  in  a  region  containing  £  =  1  and  the  curve  7  =  0  (when  it  is  to  the  left 
of  £  =  1).  Then  £0=  1  -  e(l/ct)  is  upstream  from  7  =  0,  and  i?^(£0)  +  r?T(£0)  +  co}T7.(£0) 

+  0(e)  =  0  by  (1.9).  However,  as  seen  in  section  II,  the  solution  upstream  of  the  barrier  is  of 
the  form  AFj  +  BF2  (equation  (2.9))  which  must  satisfy  (2.14)  up  to  order  e  at 
Zj  =  b(£o  -  1  -  eq).  We  then  find  (equations  (2.16)  -  (2.23))  that  IBF2I/IAF  j  |  =  0(e)  in 
[0,£q]  ,  i.e.,  outside  a  region  of  order  e  of  the  right  hand  boundary.  (The  crucial  issues  in 
the  above  handwaving  are,  of  course,  the  Taylor  series  and  the  magnitude  of  the  coefficient 
of  e  in  0(e).  If  we  try  to  consider  the  effect  of  e  0,  we  face  the  difficult  question  of  the 
dependency  of  the  series  on  e  through  equation  (3.2)!) 

CONCLUSIONS 

It  would  appear  that,  as  was  the  case  with  cj  >  Vi,  the  steady  state  solution  for  small 
flexural  rigidity  is  essentially  determined  by  the  single  boundary  condition  r?(0,r).  This  is 
certainly  true  for  equation  (3. 1 )  if  it  is  stable;  however,  neither  the  validity  of  (3. 1 )  as  a 
physical  model  nor  its  exact  relation  to  equation  (1 .2a)  has  been  rigorously  established.  The 
mathematical  difficulties  that  one  encounters  stem  from  the  presence  of  a  singular  point  in 
the  reduced  equation  and  thus  are  not  alleviated  by  a  perturbation  analysis  of  the  fourth 
order  equation. 

In  general  the  following  open  questions  remain: 

( 1 )  What  are  the  stability  properties  of  the  various  equations  studied  in  this  report? 
Are  the  solutions  stable  with  respect  to  small  changes  in  the  coefficients  (e.g.,  c^)?  This  also 
leads  to  questions  concerning  the  appropriate  class  of  functions  for  initial  conditions  and  the 
validity  of  separation  of  variables. 

(2)  What  is  the  significance,  if  any,  of  the  elliptic  region?  Does  it  simply  represent 
a  failure  of  the  model  due  to  the  approximation  q  =  (2cJ  -  1  )/2ct? 
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(3)  From  a  purely  mathematical  viewpoint,  what  is  the  relationship  between 
equations  (1.2)  and  (3.1)  for  c't<  xh\  can  one  perform  a  turning  point  singular  perturbation 
analysis,  and  does  the  limiting  solution  approach  the  solution  of  (3. 1 )? 

Finally,  we  would  like  to  suggest  that  a  reformulation  of  the  problem  with  the 
arclength  along  the  cylinder  replacing  x  as  an  independent  variable  will  probably  lead  to  a 
linear  problem  and  at  the  same  time  alleviate  many  of  the  above  difficulties  since  it  elimi¬ 
nates  the  small  angle  approximation.  Note  that  it  would  also  remove  the  inconsistency  of 
placing  the  cylinder  boundary  at  x  =  L  whereas  the  actual  boundary  must  be  at  a  point 
x  <  L  dependent  on  the  shape  of  the  cylinder  at  the  particular  point  in  time.  Such  an 
inconsistency  could  easily  be  the  cause  of  a  spurious  elliptical  region. 
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APPENDIX  A 

SEPARATION  OF  VARIABLES:  REDUCED  EQUATION 

Substituting  77  =  v(£)eRor  into  equation  (1.2a)  we  obtain 
e3hv(4)  +  b(£  -  1  -  eq)v^2)  +  (b  +  d  +  2/Jeiw)v^  ^  +  (iwd  -  eco2)v  =  0  . 
The  reduced  equation  (h  =  0)  may  be  solved  by  the  substitution  ([2] ) 


z2  =  b(£  —  1  -  eq) 


(A. 2) 


which  yields 


v77  +  —  (-1  +  —  (b  +  d  +  2/3eico))  +  v 


b  b 


fe) 


(itod  -  eu>2)  =  0  . 


,2  -  / id  w\ 

Jo  b  [b  ~e  b)  ■ 


Then  (A.3)  has  the  two  linearly  independent  solutions  ([8] ,  pg.  97) 


F1(b^z2)szl'J_i;(boz) 
F2(boZ2)  =  z"J  p(b0z) 


(A.4) 


(A.5) 


(A.6) 


where 


J„(z)  ■  ^ 

m=0 


m/  1  \  i>+2m 

(-Dm(yz) 

m!r(v  +  m  +  1) 


is  a  Bessel  function  of  the  first  kind.  (We  have  assumed  for  simplicity  that  v  is  not  an  integer; 
otherwise  we  must  also  introduce  the  Bessel  functions  Y,,  of  the  second  kind.) 

The  functions  Fj  and  F2  thus  have  power  series 


<-1>"1(P) m 

niTf-v  +  m  +  1) 


F,<z2)=  ^ 
m=0 


F  (z2j  .  f  (-')m(7z2)m(?  'V 

2  m!r(»>  +  m  +  1) 

m=0 


(A.9) 


(A.  10) 
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Their  first  and  second  derivatives  with  respect  to  are  given  by 


F*i<z2)"  2  (m- l)!r(->’+m  +  1) 


P'i<z2>=  2 

m=2 

F-2(Z2)=  £ 


F'2(z2)=  2 
m=0 


<-'>mi^(jz2P2 

(m  -  2)!  r(-i>  +  m  +  l) 

<-'>m>pjz2)mt'"' 

m!  r(v  +  m  +  1) 

ml  /i  ->\  m+p-2 

(-1)  (m  +  i>)(m  +  v -1)(-^-  z2j 


m!  r(v  +  m  +  1) 


THE  EQUATION  FOR  R0e“8^ 


For  simplicity  of  computation,  we  make  the  change  of  variable 
x  =  £  -  1  -  eq 
and  write 

v(x)  =  R(x)e"S(x)/**  , 

Substitution  of  (A.  14)  into  (2.3)  requires  the  derivatives  of  v(x)  which  are  given  by 
v'  =  [R'-g'R/Mle-S^ 


=  Ir>  .  .  £R  +  i2s.l  e-! 


,(3)  _  R(3)_  3gV  _  Vr;  _  g^R 


+  3<g)2R>  +  jgVR  (gf)3R  e-g/M 

n 2  JX2  JX3 

(4)  _  L(4)  _  4g'R<3)  _  6g”R"  __  4g(3)R'  _  g(4)R 


6(g’)2R"  12gV’Rf  3(g")2R  4g'g(3>R 


-  4(8>)3R>  _  6(g,)2g"R  +  (g>)4R  e-g/M 


(A.  18) 


Substituting  these  expressions  into  (2.3),  we  obtain  terms  in  M  and  p2. 

The  coefficient  of  p-2  is  then  set  equal  to  zero;  namely, 


(g')4Rje"8/M  +  a|X(g')2Rje~8/^  =  o 


(A. 19) 


g  =±v^x" 


which  has  the  solution 


g  =  ±  x/^T  (_x)3  ^  +  constant  . 


(A.  20) 


(A.21) 


The  function  Rg(x)  is  determined  by  equating  the  coefficient  of  the  next  highest 
term,  j pT^eT^,  to  zero: 


-4(g')3  Rq  -  6(g')2  g"R0  -  2a1xg'RQ  -  a,  xg"R0  -  a2g'R0  =  0 
Substituting  (A.20)  into  (A.22),  we  find 

-2(-x)3/2  r^|  -  ^  yrr  Rq  .  o 

Finally,  the  integral  of  (A. 23)  is  (since  x  <  0) 


l(°i  -  i\ 

RA=r„(-x)2\“>  2 / 


R0  =  ro(~x)  x  1  ' 

The  equation  for  Lq  is  the  same  as  that  for  Rq;  thus, 


A 

L0  =  ®0^~x>  2  '  1  ■ 


(A.22) 


(A.  23) 


(A.  24) 


(A.25) 
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